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ABSTRACT 

NGC 2915 is a nearby blue compact dwarf with the Hi properties of a late-type 
spiraL Its large, rotating Hi disk (extending out to i? ~ 22 _B-band scale lengths) 
and apparent lack of stars in the outer H i disk make it a useful candidate for dark 
matter studies. New Hi synthesis observations of NGC 2915 have been obtained using 
the Australian Telescope Compact Array. These data are combined with high-quality 
3.6 |J.m imaging from the Spitzer Infrared Nearby Galaxies Survey. The central regions 
of the H I disk are shown to consist of two distinct H i concentrations with significantly 
non-Gaussian line profiles. We fit a tilted ring model to the Hi velocity field to derive a 
rotation curve. This is used as input for mass models that determine the contributions 
from the stellar and gas disks as well as the dark matter halo. The galaxy is dark- 
matter-dominated at nearly all radii. At the last measured point of the rotation curve, 
the total mass to blue light ratio is Mioi/Lb ~ 140 Mq/Lb, making NGC 2915 one 
of the darkest galaxies known. We show that the stellar disk cannot account for the 
steeply-rising portion of the observed rotation curve. The best-fitting dark matter halo 
is a pseudo-isothermal sphere with a core density po ^ 0.17 ± 0.03 Mq pc~^ and a 
core radius Tc ^ 0.9 ± 0.1 kpc. 
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INTRODUCTION 



■ Dwarf galaxies are thought to d ominate the cosmic 
\ scenery in terms of number de nsity l|Springel et al.l l2005l : 
! ISandage et afll 19851 : lMatedll998l ') and therefore form an im- 
' portant morphological class of galaxies. Several investiga- 
tors have studied the mass dist r ibutions of these galajcies 
(e.g. [Carignan fc Beaulieul 19891: Cote et al. 1991 ;_ Broeil j 



Il992l: iMeurer et al.l 1199^ Ide Blok fc McGaughl Il997l. and 



references therein). ISwatersI (199^) was the first to study 
the dark matter (DM) properties of a large representative 
sample of nearby dwarf galaxies as part of the Wester- 
bork Hi Survey of Spiral and Irregular Galax;ies (WHISP, 
ISwaters et~all 120021 ). More recently, a morphologically di- 
verse sample of nearby galaxies, including dwarf systems, 
was observe d as part of The H l Nearby Galaxy Survey 
(THINGS, I Walter et all |2008| ). These cumulative efforts 
have led to an understanding that DM plays an important 
role in the d ynamics and evolution of dwarf galaxies. 

ISwatersI 11999) found that none of the observed rota- 
tion curves in his sample of late-type dwarf galaxies decline 
at outer radii. This was further confirmation that these sys- 
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tems have a large dark mass fraction. Dwarf galaxies there- 
fore serve as ideal candidates to test theories of DM. Ob- 
servational studies of the mass distributions of these dwarf 
systems have shown their DM halos to be well-modelled 
by a pseudo iso-ther mal sphere with an approximately 
constant-density co re (|de Blok et al ] | 19961 : ICote et al. I I2OOOI : 
Ide Blok et al.lfioO^ '). Theoretically, the study of the prop- 
erties of DM halos in a hierarchical clustering context is 
carried out mainly in the f orm of Cold Dark Ma tter (CDM) 
numerical simulations (e.g . iNavarro et al.„1996, : , Moore et al.l 
ll999l : ISpringel et al.|[200a '). The equilibrium density profile, 
p(r), of simulated DM halos is found to vary wit h radius, 
r, as p(r) cx r" with a ~ —1 (|Navarro et al.|[l997l h thereby 
predicting extremely steep inner density profiles. This dis- 
crepancy between the observed shapes of DM halos and the 
predictions from numerical simulations has become known 
as the "cusp/core" problem. 

Attempts at reconciling observations with theoretical 
predictions depend crucially on accurate dynamical analy- 
ses of nearby galaxies. To test DM halo models one wants 
galaxies that are as dark-matter-dominated as possible. In 
this respect, a potentially useful candidate for DM stud- 
ies is NGC 2915. This galaxy is classified as a nearby 
{D ~ 3.78 Mpc, iKarachentsev et al.l 120041 ') blue com- 
pact dwarf according to its optical appearance, yet has the 
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Hi morpholo g y of a late-type spiral. Detailed imaging by 
[m eurer et al. I |l994l ) showed that the optical appearance of 
NGC 2915 is dominated by two main stellar populations: 
a compact blue population, which is the location of on- 
going high-mass star formation; and a more diffuse, older 
red population. What makes this galaxy so interesting is 
the fact that, when observed at 21 cm, the stellar core of 
NGC 2915 is seen to be completely embedded in a huge 
Hi disk extending out to i? ~ 22 _B-band scale lengths 
l|Meurer et al.|[r996l ). Furthermore, this gas disk has well- 
defined spiral structure that is not seen in the stellar disk. 
No significant star-formation is observed in the outer parts 
of this gas disk. Only a few faint HII regions were detected 
bv lMeurer et al.l (|l999t ) after carrying out deep Ha imaging 
of the disk. 

The Hi disk of NGC 2915 serves as an ideal tracer of 
the gravitational potential out to radii of _R ~ 10 kpc, far 
beyond the visible radial extent of the stellar disk. Due to 
the apparent lack of stars in its outer parts, the dynam- 
ics of the H I disk are an a lmost direct tracer of the DM 
distribution. iMeurer et al] (|l996l ') targeted NGC 2915 for 
DM studies and found the galaxy to be DM-dominated 
at nearly all radii with a total mass to B-band light 
ratio Mtot/Ls -c^ 76 Mq/Ijq (assuming a distance of 
5.1 Mpc). They also found the galaxy to have a dense 
(po ~ 0.1 Mq pc~"^) and compact (rc ~ 1 kpc) DM core. 

Several investigators have attempt ed to explain t he ob - 
served H i distribution of NGC 2915. [Bureau et al] (|l999l ') 
studied the dynamics of the central Hi region and the spi- 
ral pattern of the outer disk. For both they calculated a 
common, slow pattern speed of Q,p = 8.0 ± 2.4 km s~^ 
which they associated with the figure rotation of a tri- 
axial DM halo. They also proposed that some DM is dis- 
tributed in the disk of NGC 2915, thereby making it grav- 
itationally u nstable to the f ormation of the observed spi- 
ral structure. iMasset &: Bureaul (j 20031 ). using hydrodynam- 
ical si mulations, further explored the ideas of [Bureau et al.l 
l|l999l ). They showed that the observed spiral structure can 
be accounted for by either an unseen bar or a rotating 
tri-axial DM halo. However, the mass of the required bar, 
Mi,ar ~ 5 X 10^ Afo, is very large in comparison to the 
total stellar mass, thereby making the nature of such a bar 
problematic. The required pattern speed of the tri-ax;ial DM 
halo is s ignificantly larger than those from numerical sim- 
ulations. iMasset fc Bureaaa (|2003i l disfavoured the external 
perturber scenario. They also found that while a heavy disk 
is able to account for the main features of the observed H I 
morphology, it fails to match the observed gas dynamics. A 
satisfactory explanation for the various morphological and 
kinematic features of the Hi in NGC 2915 is therefore still 
lacking. 

In this paper we investigate the Hi dynamics and DM 
distribution of NGC 2915. We use data from new Hi syn- 
thesis observations of NGC 2915 carried out using the Aus- 
tralian Telescope Compact Array as part of the Southern 
Hemisphere extension of The H I Nearby Galaxy Survey. Our 
H I observations are significantly deeper and have better spa- 
tial resolution than any other H I observations of NGC 2915. 
These data are complemented by high-quality 3.6 \im in- 
frared observations of the stellar disk, carried out as part 
of the Spitzer Infrare d Nearby Galaxies Survey (SINGS, 
iKennicutt et al.ll2003l ). The Hi observations are presented 



in Sec. [2] while the Hi data products appear in Sec.|3] This 
paper focuses on the regular gas dynamics of NGC 2915. In 
Sec. m a rotation curve out to i? ~ 9.3 kpc is determined. 
The far-infrared observations of the stellar disk are com- 
bined with our new Hi observations in Sec. [5] to produce 
a mass model that determines the contributions at various 
radii of the stars, gas and DM. We use two different pa- 
rameterisations of the DM halo to reconstruc t the observed 
rotat ion curve, one of which is the NFW halo (iNavarro et al.l 
\l99i ) favoured by numerical simulations while the other is 
the observationally motivated pseudo-isothermal sphere. Fi- 
nally, in Sec. [S] we summarise our results and present our 
conclusions. 



2 H I OBSERVATIONS AND DATA 
REDUCTION 

2.1 Data acquisition 

NGC 2915 was observed between 23 October 2006 and 
2 June 2007 (project number C 1629) with six different 
ATCA configurations using all six antennas. A single run 
consisted of a primary calibrator observation, regular sec- 
ondary calibrator observations and source observations. The 
calibrator sources PKS 1934-63 and PKS 0823-500 were 
used as primary and secondary calibrators respectively. Each 
of the EW-352, 750D, 1.5B and 1.5C runs was approximately 
12 hours long while 24 hours were spent in the 6A configu- 
ration. Besides these data, archival data were also incorpo- 
ratecQ. iMeurer et al] l|l996l ) determined the Hi diameter of 
NGC 2915 to be ~ 0.32°. The Hi extent of NGC 2915 there- 
fore falls well within the ~ 0.54° field-of-view of the ATCA 
dishes when observing at 21 cm. For our observations, the 
telescope pointing centres were set to the optical centre of 
NGC 2915 (a2ooo = 09'^ 26"" 11.5% c^zooo = -76° 37' 35") with 
no mosaicking required. The correlator was set to use 512 
channels with a bandwidth of 8 MHz, centered at 1418 MHz. 
The resulting velocity range is —335 km s~^ to 1353 km s~^ 
with an approximate channel spacing of 3.2 km s~^. Table 
[1] provides a summary of all our observing setups. 

2.2 Hi Data cubes 



The MIRIAD software package (jSault et al.lll995l ) was used 
to take the raw uv data from the correlator through to the 
image analysis stage. All velocities are measured relative to 
the barycentre rest frame. The first five and last five corre- 
lator channels were flagged. The data were then split into 
primary calibrator, secondary calibrator, and source subsets. 
The required corrections to the antenna gains, delay terms 
and bandpass shapes were determined. The time-varying 
phases and antenna gains were calibrated based on obser- 
vations of the secondary calibrator. 

The calibrated source data were continuum-subtracted 
by fitting and subtracting a flrst-order polynomial to the 
line- free channels which, for our observations, were chan- 
nels 50 — 100 and 350 — 460. Having isolated the Hi signal, 
image cubes were produced from the continuum-subtracted 
uv data. Data cubes were produced using natural weighting 

^ Project number C191, principle investigator: Meurer. 
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Table 1. Summary of NGC 2915 observing setups 



1 


2 


3 


4 


5 


Conf. 


Date 


Start 


End 


Dur. 




(yy-mm-dd) 


(hh-mm-ss) 


(hh-mm-ss) 


(hrs) 


EW352 


2006-10-23 


11:28:55 


23:17:45 


10.39 


750D 


2007-03-14 


04:58:45 


15:56:10 


9.64 


1.5B 


2006-11-24 


19:29:05 


06:51:35 


9.95 


1.5C 


2007-05-02 


04:37:25 


15:36:05 


9.64 


6A 


2007-02-10 


10:34:35 


21:32:25 


9.64 


6A 


2007-02-16 


05:12:35 


00:42:15 


17.06 



Comments on columns: Column 1: ATCA configuration used; 
Column 2: date of observation (UT); Column 3/4: start/end of 
observations (UT); Column 5: time on source. 

as well as robust weighting with a Brigg's visibility weight- 
ing robustness parameter of 0.2. Each of the d irty images 
was d econvolved using a Steer Clean algorithm (|Steer et al.l 
Il984h to produce an output map of clean components. Chan- 
nels were either cleaned down to a flux cut-off of 2.5 times 
the typical r.m.s. of the flux in a line-free channel or for 50 
000 iterations, whichever condition was met first. The typ- 
ical number of iterations was ~ 20 000 - 30 000. After the 
deconvolution process, each of the clean components was 
convolved with a Gaussian approximation of the dirty beam. 
The beam's full width at half maximum was 17" x 18.2" and 
10.2" X 10.2" for the naturally-weighted (NA) and robust- 
weighted (RW) data cubes respectively while the channel 
width was set to dV = 3.49 km s~^ for both. No Hanning 
smoothing was applied to either data cube. For each cube, 
the noise in a line-free channel is Gaussian distributed with 
a mean of /i ~ —0.02 mjy beam~^ and a standard devi- 
ation of (J ~ 0.60 mJy beam~^ for the NA cube, while 
fj, ~ —0.004 mJy beam~^ and a ~ 0.58 mJy beam~^ for 
the RW cube. 



2.3 Hi Moment maps 

All Hi moment maps were generated using the Gronin- 
gen I mage Processing System (GIPSY, Ivan der Hulst et al.l 
Il992l ') together with the NA Hi data cube. The entire Hi 
data cube was smoothed, using a Gaussian convolving func- 
tion, to a resolution of 34" x 36.4". A flux cut-off was ap- 
plied to the smoothed cube at 2.5a. For each channel of the 
smoothed data cube, any flux that was clearly not spatially 
connected to, or associated with the galaxy emission was 
flagged and removed. The blotted smoothed cube was then 
applied as a mask to the original 17" x 18.2" cube. 

The Hi total intensity and 2nd-order moment maps 
were constructed from the zeroth- and second-order mo- 
ments of the data cube respectively. A map in which each 
pixel is equal to oy/N , where N is the number of un-blanked 
channels for a given line profile that would contribute to the 
moment maps, was created. Such a construction was per- 
mitted since no Hanning smoothing was applied to the H I 



data cube. Dividing each pixel in the total intensity map by 
ay/N produced a dimensionless signal-to- noise (S/N) map. 
The average intensity of all the pixels in the total intensity 
map that had a corresponding signal-to-noise ratio in the 
range 2.75 ^5 S/N 5C 3.25 was used as an intensity cut-off for 
the total intensity map. This removed, in a robust manner, 
pixels from the outer disk that had low S/N values. This new 
H I total intensity map was then used as a mask on other H I 
maps. 

The method of using Gauss-Hermite polynomials to 
parameterise line profi l es of galaxies was introduced by 
Ivan der Marel fc Franxl (|l993l '). While Gaussians are rea- 
sonable approximations to many realistic line profiles, 
higher order Gauss-Hermite polynomials are able to better 
capture profile asymmetries and non-Gaussian deviations. 
The Hermite method has already been successfully imple- 
mented by various auth ors (van der Marcl_& Franx 199^; 
iNoordermeer et all [2OO7I : Ide Blok et al.ll200d '). We fitted a 
third-order Gauss-Hermite polynomial to each NA H I data 
cube line profile to generate a so-called Gauss-Hermite H I 
velocity field. Three filters were used simultaneously when 
fitting the profiles: 1) profiles with fitted peak fluxes below 
2.5a were excluded; 2) proflles with a fitted line width less 
than the channel width were excluded; and 3) fitted profile 
peaks had to be within the velocity range of the data cube. 



3 HI PROPERTIES 
3.1 Channel maps 

A grey-scale representation of the 17" x 18.2" resolution 
channel maps is shown in Fig. [1] cr ~ 0.6 mJy in all chan- 
nels. H I structures are visible at large and small scales. At 
large scales, the emission clearly exhibits the usual pattern of 
a rotating disk. The Hi disk of NGC 2915 exhibits clear spi- 
ral structure (c.f. Hi total intensity map. Fig. On small 
scales, the gas in the central regions forms two dominant 
over-densities in nearly all of the channels, with fainter emis- 
sion in-between them. 



3.2 Global profile 

A global H I profile, extracted from the blotted NA H I data 
cube by summing the H I emission in each velocity channel, 
is presented in Fig.O Using a synthesis array to estimate the 
total H I flux of a galaxy can result in an underestimate due 
to missing short baselines. The shortest baseline used for our 
observations was 31 m (as part of the EW352 antenna con- 
figuration) , implying that the ob servations are less se nsitive 
to structures larger than ~ 23' . iMeurer et al.l (|l996l ). how- 
ever, estimated the H I extent of NGC 2915 to be ~ 19' . The 
amount of H I missed in the global H I profile is therefore ex- 
pected to be small. In Fig. [5] we compare our determination 
of the global H i profile for NGC 2 915 t o the HIPASS global 
Hi profile from iKoribalski et all (|2004l ). The HIPASS data 
are calibrated with 8' spatial binning and 13 km s~^ chan- 
nel width. Due to the large HIPASS beam, as well as the 
fact that the Parkes telescope is a single dish telescope, the 
amount of flux missed in the HIPASS observations should be 
negligible. It is clear from Fig. [2]that no significant amount 
of flux has been missed by our observations. 
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Figure 1. Channel maps of NA Hi data cube. Every second channel is shown. The heliocentric radial velocity (km s~^) of ea<;h channel 
is shown in the upper left corner. The half-power-beam-width is shown in the bottom left corner. Grey-scale range is from -0.6 mjy 
beam~^ to 10.0 mJy beam~^. The r.m.s. noise in a channel is ~ 0.6 mJy beam"^. 



The resulting global Hi profile is double-horned and 
asymmetric. EYom this profile we derived the profile widths, 
the systemic velocity, and the total Hi mass. The profile 
widths were calculated as the differences between the high 
and low velocities of the galaxy at 20% and 50% of the 
peak flux density {Vf;^l, V?°I°, V^°^^, Vf°I° respectively). 
The profile widths were not corrected for inclination. Two 



systemic velocity estimates were determined: 1) from the 
global profile, at the midpoint of the two velocities at the 
50% level, 2) from the equation 

Vsys = 0.25 X {V^t + + + (1) 
Finally, the total H I mass was calculated as 
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Figure 1. Continued. 



Mhi = 2.36 X 10^ X X 



FdV, 



(2) 



where D is the distance to the galaxy in Mpc and J FdV 
is the total Hi line flux in Jy km s~^. This determination 
of the total Hi mass assumes that the Hi is optically thin. 
Throughout this paper we adopt the D ~ 3.78 Mpc dis- 
tance determination of iKarachentsev at al.l ()2004l l. Table [2] 
lists our determinations of the above mentioned quantities . 
Assuming a distance of D = 5.1 Mpc, iMeurer et al.l (| 19961 ) 



measured a total Hi mass of AIhi ~ 1.27 x 10 AIq for 
NGC 2915. Using a distance of D = 3.78 Mpc, their total 
mass estimate reduces to Mhi ~ 7 x 10* Mq, consistent 
with our estimate of Mhi ~ 4.4 x 10* Mq. 



3.3 Total intensity map 

The full H I total intensity map extracted from the NA H I 
data cube is shown in the upper panel of Fig. |3l The edge of 
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Table 2. Quantities derived from the global Hi profile. 



1 

Mhi 
(10" Mq) 


2 

W50 
(km s~l) 


3 

(km a~l) 


4 

(km s~l) 


5 

T^Eqn.l 

"ays 
(km s~l) 


4.4 


186.8 


170.0 


460.5 


460.8 



Comments on columns: Column 1: total Hi mass; Column 2: 
velocity width at 50% of the peak flux; Column 3: velocity width 
at 20% of the peak flux; Column 4; systemic velocity from W50 
midpoint; Column 5: systemic velocity from Egn. fTI 
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Figure 2. Global Hi profile from our observations (solid curve) 
and HIPASS (dashed curve). The dotted vertical line shows the 
systemic velocity as calculated using Eqn[T] 

the Hi disk (the outer-most contour level) lies at J? ~ 510" 
and therefore extends out to i? ~ 20 B-band scale lengths. 
For comparative purposes the IRAC 3.6 \im image of the 
stellar disk is shown in Fig. |4l To determine the photomet- 
ric centre we fitted ellipses to three 3.6 vim flux density an- 
nuli near the edge of the old stellar disk. A 3.6 [im flux 
density contour at a level of 1.2 MJy ster~^, the average 
flux density of the second annulus, is shown in Fig. |4l The 
average flux densities of the other two annuli are 0.8 and 
1.7 MJy ster"^ The average centre position of the ellipses 
fitted to the fiux density annuh, 02000 = 09*" 26™ 12.611% 
fcooo = —76° 37' 37.80", was used to estimate the position 
of the photometric centre (shown as a cross in Fig. |4]). For 
comparison, the 2MASS determination of the centre of the 
system is 02000 = 09*' 26™ 11.53", (52ooo = -76° 37' 34.80" 
l|Skrutskie et al.l |2006| ). To facilitate length-scale compar- 
isons between IR and Hi images, an ellipse with semi-major 
axis length a = 50' ', roughly equal to t he 7i-band R25 ra- 
dius determined bv lMeurer et al.l (|l994l ). is included in the 
figures. In thei r study of the Hi c ontent of 73 late-type 
dwarf galaxies, ISwaters et al.l (|2002l ) found that the ratio 
of the Hi extent to the optical diameter, defined as 6.4 disk 
scale lengths, is 1.8 ± 0.8 on average. The same ratio for 
NGC 2915 is ~ 3.4, close to double that value. 

Our high-resolution H I observations resolve the central 
regions of NGC 2915 into two distinct Hi concentrations 
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Figure 3. Grey-scale plot of the Hi total-intensity map (top 
panel) with a zoom-in of the central region (bottom panel). In 
both panels the intensity scale runs from 5 mjy beam~^ (white) 
to 0.2 Jy beam~^ (black). Contour levels run from 10 - 
60 mJy beam~^ in steps of 10 mJy beam~^ and from 80 - 
230 mJy beam~^ in steps of 20 mJy beam"'^. The position of the 
photometric centre estimated using the 3.6 ]im image is marked 
with a cross. The hatched circles in the lower corners represent 
the half power beam width of the synthesised beam. The black 
and white ellipses have a semi- major axis length a = 50", an incli- 
nation i = 55°, a position angle PA = 317° and a centre position 
equal to our determination of the photometric centre. The solid 
black lines represent the position- velocity slices shown in Fig. [9] 

(Fig. [31 lower panel). Their centres are separated by ~ 60" 
(1.1 kpc) and have a combined mass of ~ 1.9 X 10'' Mq. Also 
visible is a plume-like H l feature to the North- West of these 
central Hi concentrations with a mass of ~ 2.8 X lO'^ Mq. 
The observed spiral structure of the outer disk is asym- 
metric: the arm beginning at the North- Western central H I 
concentration can be traced over ~ 200° in azimuthal an- 
gle while the arm beginning at the South-Eastern Hi con- 
centration can be traced for only ~ 120° . From the H I 
total intensity map an azimuthally averaged, inclination- 
corrected Hi surface density profile was constructed. Each 
ring has a width of dR = 17", a position angle (measured 
anti-clockwise from North to the receding major axis) of 
PA = 285° and an inclination of i = 55°. The profile is 
shown in Fig. [5] 

In an attempt to better constrain the Hi distribution 
of NGC 2915, an Hi total intensity map was extracted from 
the RW H I data cube. For the sake of brevity this map is not 
shown here, yet it suffices to say that the higher spatial res- 
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Figure 4. Grey scale plot of the 3.6 |J.m IRAC Spitzer im- 
age. The intensity scale runs from -0.2 MJy ster~^ (white) to 
2.0 MJy ster~^ (black). The single white contour is at a level of 
1.2 MJy ster~^ and represents the average flux of one of three 
flux annuli used to estimate the position of the photometric cen- 
tre. The estimated photometric centre is marked with a cross. 
The black ellipse is the same as that shown in Fig. |3] 
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Figure 5. Hi surface density radial profile. All surface densities 
have been inclination-corrected. Error bars represent the r.m.s. 
spread of the surface densities in each azimuthally averaged ring. 



olution confirms the above-mentioned nature of the central 
Hi concentrations. However, the very low S/N ratio in the 
outer parts of this map yields it unsuitable for the study of 
low column density H I features, with only the central- most 
Hi features standing out about the noise. This map does 
not allow for a more accurate study of the H I distribution 
and is therefore not discussed further. 



3.4 Velocity field 

The third-order Gauss-Hermite velocity field extracted from 
the NA H I data cube is shown in Fig. (6) The blank pixels in 
this map are those that were rejected by the fitting routine 



according to the criteria mentioned in Sec. l2.3llde Blok et al.l 
(2008) showed that the Gauss-Hermite parameterisation of 
a line profile gives a robust estimate of the peak-associated 
velocity of the profile. Inner contours running approximately 
parallel to the kinematic minor axis are caused by V(ii) oc R 
rotation while the outer contours running radially away from 
the centre are due to a constant rotation speed. The overall 
"S-shaped" distortion of the contours is indicative of a kine- 
matic warp within the disk. Small wiggles along the outer 
contours are caused by streaming gas motions. The sharp 
kinks at inner radii suggest the presence of non-circular ve- 
locity components within the gas. Finally, differences in the 
shapes of the iso-velocity contours on the receding and ap- 
proaching halves of the galaxy suggest a certain degree of 
kinematic lopsidedness within the galaxy. 

A third-order Gauss-Hermite velocity field was also ex- 
tracted from the RW Hi data cube. Due to the high noise 
levels in the RW Hi cube, most line profiles were rejected 
by the fitting filters and hence not fitted. The general char- 
acteristics of the fitted peaks of the remaining profiles are 
very similar to those of the NA velocity field, differing on 
average by ~ 5 km s~^. Rather than showing this Hi veloc- 
ity field here, a tilted ring model is fitted to it in Sec. l4.3l to 
demonstrate that the resulting rotation curve is very similar 
to the one derived using the NA third-order Gauss-Hermite 
velocity field. Motivated by this result is the decision to use 
only the NA velocity field for our kinematic analyses. 



3.5 Second-order moment map 

The second-order Hi moment map extracted from the NA 
H I data cube is shown in Fig. [T] The highest second-order 
moments are found near the stellar core of the galaxy. The 
radial profile of the second-order moment map is shown in 
Fig. [3 K the second-order moments are interpreted as a mea- 
sure of the gas velocity dispersion, then the average outer 
disk dispersion of agas ~ 10 km s~^ is close to the constant 
gas velocity disper sion value of a gas = 11 km s~^ found by 
iLeroY et al] l|2008l ) for a sample of THINGS dwarf galaxies. 



3.6 Non-Gaussian line profiles 

Position-velocity slices taken though the centre of the NA 
Hi data cube betray the presence of broad and multi- 
component line profiles at inner and outer radii. Figure [9] 
shows two position-velocity slices, the positions of which are 
shown in the Hi total intensity map (Fig.O top panel). The 
first slice was taken at a position angle of 322° , through both 
of the central H I concentrations. This slice suggests that the 
galaxy contains a fast-rotating gas component near its cen- 
tre, clearly manifesting itself as a sharp velocity spike with a 
peak velocity of 1/ ~ 550 km s~^. Line profiles near the cen- 
tre of this position-velocity slice show significant departures 
from Gaussianity. 

This position-velocity slice also reveals asymmetric line 
profiles on the receding side of the galaxy, seen in the upper- 
left quadrant {V > Vsys ~ 465 km s~^ , angular offset < 0) as 
an "H I beard" that is lagging in velocity relative to the main 
disk. The beard emission spatially coincides with the plume- 
like Hi feature seen in the Hi total intensity map. The esti- 
mated mass of this morphological feature, ~ 2.8 x 10^ Mq, 
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Figure 6. Top panel: Grey-scale plot of the third-order Gauss- 
Hermite H I velocity field. Bottom panel: A zoom-in of the central 
region. The intensity scale runs from 300 - 600 km s~^. Contours 
are separated by 15 km s~^ with the thick contour marking the 
systemic velocity at 465 km s~^ derived by fitting a tilted ring 
model to the velocity field. White and black contours represent 
receding and approaching halves of the galaxy respectively. The 
hatched circles in the lower corners represent the half power beam 
width of the synthesised beam. The ellipse is the same as that 
shown in Fig. l4l 
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Figure 7. Grey-scale plot of the second-order H I moment map. 
The grey scale runs from 6 km s~^ (white) to 30 km s~^ (black). 
Contour levels are at 6, 9, 12, 18, 22, 26, 28, 30 km s"!. The 
hatched circle in the lower left corner represents the half power 
beam width of the synthesised beam. The ellipse is the same as 
that shown in Fig. |4] 
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Figure 8. Second-order Hi moment map radial profile. Error 
bars represent the r.m.s. spread of the 2nd-order moments in each 
azimuthally averaged ring. 



is 5.6% of total H I mass. One interpretation of this kinemat- 
ically anomalous gas component is that it is associated with 
pristine gas being accreted from the nearby inter-galactic 
space onto the outer disk of NGC 2915. If confirmed, this 
would be the first evidence of such an event in a dwarf sys- 
tem like NGC 2915. 

The second position-velocity slice was extracted almost 
parallel to the minor axis, through the South Eastern central 
Hi concentration which contains many double-component 
profiles. Line profiles in this central region are split by 
~ 30 km s~^ on average, implying that the sharp rise 
seen toward the centre of the second-order H I moment map 
(Fig. [7]) over-estimates the true H I velocity dispersions of the 
inner disk. Since the disturbed Hi line profiles are located 
within R ~ 150" of the centre of NGC 2915, it is plausible 
that the central gas dynamics of NGC 2915 are largely dic- 
tated by the stellar winds of its young stellar population. In 
a forthcoming paper (Elson et al. 2010, in prep.) we present 
detailed modeling results of the central H I dynamics. 



4 H I DYNAMICS 

The H I line data allow us to study the H I gas dynamics. In 
this section we derive a rotation curve for NGC 2915. 



4.1 Tilted ring model 

The standard method of producing a rotati on curve involve s 
fitting a tilted ring model to a velocity field (|Begemanlll987l ) . 
NGC 2915's intermediate inclination of i ~55° makes it a 
well-suited candidate for such an analysis. The method in- 
volves modeling the galaxy as a set of concentric rings within 
which the gas orbits about the kinematic centre. Each ring 
has a set of defining parameters: central coordinates and 
Yc, inclination i, systemic velocity Vsya, position angle PA 



The dark matter content of the blue compact dwarf NGC 2915 9 



500 



E 



o 
> 



400 



500 



E 



o 
_o 

CD 
> 



400 




PA in the sky plane by 



+ 0.1 -0.1 

Angular offset (degrees) 

Figure 9. Position-velocity diagrams extracted along 1) a po- 
sition angle of 322° (top) and 2) a position angle of 220° (bot- 
tom) . The slice positions are shown in the Hi total intensity map 
(Fig. [Sj. The thickness of the slices is that of a single pixel in the 
data cube (2.5") and the velocity resolution is 3.5 km . The 
intensity scale runs from — Scr to (white) 40o- (black). Contours 
start at 2cr in steps of 3cr. The dashed horizontal line marks the 
systemic velocity at Vsys = 465 km . 



(in the sky plan^EI), and the rotation velocity Vrot of the ma- 
terial within the ring. In the case that expansion velocities 
within the disk are ignored, the standard algorithm carries 
out a least squares fit to 



Vios{x,y) = Vsya + Kot sin(i) cos(6l), 



(3) 



where Vios is the line-of-sight velocity, x and y are rectangu- 
lar coordinates on the sky and 9 is the angle from the major 
axis in the galaxy plane. 9 is related to the position angle 



Measured anti-clockwise from North to the receding semi-major 
axis. 



cos(6') = 
sm(9) = 



-(a; - Xc) sin(Pyl) + (y - Yc) cos(PA) 
r 

-{x- X,) cos(PA) + {y- Yc) sin(P^) 



(4) 
(5) 



r cos(i) 

where r is the radius of the ring i n the galaxy pla ne. 

The GIPSY task ROTCUR (|Begemanlll987l ) was used 
to fit tilted ring models to the third-order Gauss-Hermite 
H I velocity field extracted from the NA H i data cube. Rings 
width of dr = 17" were used to ensure that adjacent rings 
were largely independent of one another. When fitting Eqn.[5] 
to the data, each datum was weighted by jcos(6)| so that 
points closer to the major axes held more weight. All points 
within 10° of either side of the minor axes were excluded 
from the fit. Both sides of the galaxy were used for the fit- 
ting. 



4.2 Fitting procedure 

Each tilted ring should be centred on the dynamical cen- 
tre of the galaxy, assuming it is constant, which is usu- 
ally estimated by the least-squares fitting algorithm to- 
gether with the other tilted ring parameters. In the case 
of NGC 2915, however, the kinematic centre determinations 
from ROTCUR varied with radius in a non-systematic man- 
ner. We therefore used the photometric centre of the 3.6 vim 
emission as the kinematic centre. Using the photometric and 
kinematic centres interchangeably in this manner is reason- 
able if one assumes that the stars lie at the bottom of the 
gravitational potential. We carried out a single ROTCUR 
iteration with all parameters allowed to vary freely. The re- 
sulting Xc and Yc values are shown as filled black circles 
in panels A and B respectively of Fig. 1101 These fitted Xc 
and Yc positions deviate on average by 8.3" ± 0.3" and 
7.2 " ± 0.25" (152 pc and ~ 132 pc) respectively from 
the photometric centre, thereby placing them well within 
a single half power beam w idth. This average deviat ion is 
consistent with the results of lTrachternach et al.l (|2008l ) who 
found that approximately 50% of their kinematic centres de- 
rived for ~ 1000 individual tilted rings fitted to the THINGS 
Hi velocity fields differed by less than a beam width from 
their best centre estimates. In the case of the THINGS sam- 
ple, one beam width corresponds to ~ 200 pc on average. 
Throughout the tilted ring fitting procedure, Xc and Yc were 
therefore each fixed to the position of the photometric centre 
(solid lines in panels A and B of Fig. llOp 

The remaining tilted ring parameters (Vsys, PA, i, Vrot) 
were derived iteratively. We first allowed all parameters 
to vary freely. This gave a general feel for the parame- 
ter behaviour with r. Little scatter was seen in the radial 
run of Vsys, which was then fixed to the average value of 
465 km s~^ in subsequent iterations (compare with the es- 
timate of Vsys = 460.8 km s~^ derived from the global Hi 
profile using Eqn.[l|. The PA, i, and Vrot parameters were 
allowed to vary freely in the subsequent iterations. The PA 
and i runs showed relatively little scatter except for a clear 
change from outer to inner disk beginning at 7? ~ 240" where 
the inclination value changed from i ~ 55° to i ~ 75° . Our 
kinematically derived position angles were consistent with 
the position angles of ellipses fitted to flux contours in the 
Hi total intensity map. Several more iterations with various 
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Figure 10. Radial variations of tilted ring parameters fitted to the third-order Gauss-Hermite velocity field. Panels are A: central X 
positions, B: central Y positions, C: position angles, D: inclinations, E: systemic velocities, and F: circular rotation velocities. Filled 
circles show the radial variations just before the parameter profiles were smoothed and fixed in order to derive the final rotation curves. 
The smoothed profiles used to construct models SI and CI are shown as solid curves except for the model CI inclination profile which 
appears as a dashed curve (offset by -2° for the sake of clarity). Panel F shows the final rotation curves derived for models SI and CI 
(open and filled squares respectively). For the sake of clarity the rotation curves have been offset by 8.5". Error bars represent the r.m.s. 
spread of the velocities in a given ring. 



combinations of fixed and free parameters were carried out 
in order to check the stability of the radial parameter pro- 
files. The sharp rise in i for R < 240" always remained with 
little scatter. 

For the final iteration, the Xc, Yc and Vsys parameters 
were fixed to values that were constant with r while the 
inclination and position angle radial profiles were smoothed 
and then fixed. The smoothed versions of the profiles are 
shown as solid curves in panel A, B, C, D, E in Fig. 1101 while 
the filled circles represent the parameter values just before 
they were smoothed and fixed. Smoothing the PA and i 
radial distributions did not affect the fina l rotation curve in 
any noticeable manner. As d e Blok et al.l (2008) point out, 
the smoothing only affects the point-to-point scatter of the 
profiles and in no way affects the resolution of the radial 
distributions. To derive the final rotation curve, only Vrot 
was allowed to vary. One of the main sources of uncertainty 
in the derivation of the inner rotation curve is the rise in 
inclination of ~ 20° towards inner radii preferred by the 
tilted ring models. It is not clear whether this rise is due to a 
genuine kinematic warp, or whether it is merely an artefact 
of the non-Gaussian line profiles near the galaxy's centre. 
We therefore derived two final tilted ring models, one with 
and one without the rise in inclination (solid and dashed 
curves in panel D of Fig.[lO]respectively). Both models used 
the same X^, Yc, Vsys and PA radial profiles. We hereafter 
refer to these two tilted ring models as models SI (Steep 
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Figure 11. Rotation curves from tilted ring model SI derived 
separately for the approaching (filled circles) and receding (open 
circles) sides. For the sake of clarity the rotation curve of the 
approaching side has been offset by 8.5". Error bars represent 
the r.m.s. spread of the velocities in a given ring. 



Inclination) and CI (Constant Inclination). In panel F of 
Fig. 1101 the rotation curves of both models are presented. In 
each of the panels A - E in Fig. IIOI the error bars represent 
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Figure 12. Position-velocity slices though the kinematic centre of NGC 2915 at various position angles. Overlaid on each slice is the 
rotation curve of tilted ring model SI projected using the smoothed parameter radial profiles as shown in Fig. 1101 The thickness of the 
position-velocity slices is that of a single pixel in the data cube (2.5") and the velocity resolution is 3.5 km s~^. The dashed horizontal 
line shows the systemic velocity at 465 km s~^ adopted for our kinematic analyses. Contours start at 2cr in steps of 2.5(t. 



the formal least squares errors from the ROTCUR fitting 
routine. In panel F, the error bars represent the r.m.s. spread 
of the velocities in a single ring. 

4.3 Rotation curve 

The final tilted ring models fitted to the third-order Gauss- 
Hermite velocity field are presented in Fig. 1101 While the 
stellar disk seen in the 3.6 \im image (Fig. |4]) extends as far 
as i? ~ 100", the inner portion of the rotation curve rises as 
V{R) oc R out to 7? ~ 150". Beyond this radius the rotation 
velocity remains almost constant. For model SI, the average 
rotation velocity for R ^ 187" is V = 81.9 ± 1.6 km s~\ 
Despite using different inclination profiles to derive the two 
final rotation curves, the maximum absolute difference be- 
tween them isV = 9.3±12.2 km s"^ at i? = 136". The aver- 



age absolute difference for R ^ 187" is V = 4.3±4.8km s"\ 
The resulting two rotation curves are thus very similar to 
one another. For each tilted ring model we constructed sep- 
arate rotation curves for the approaching and receding sides 
of the galaxy. In the case of model SI (Fig. Ilip . the ro- 
tation curve differs significantly between the two sides of 
the galaxy. Within the steeply rising portion, differences of 
~ 20 km s~^ are observed as wel l as at J? ~ 290". The same 
is true in the case of model CI. ISwaters et al] (|l999l ) state 
that galaxies with asymmetric global profiles often have ro- 
tation curves that are more slowly rising on one side of the 
galaxy than the other. This statement seems true in the case 
of NGC 2915. 

To check that the two tilted ring models are consis- 
tent with the data, we used the smoothed radial parameter 
profiles to project the rotation curves onto various position- 
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Figure 13. The model SI rotation curve (open circles) compared 
to the rotation curve derived by fitting a tilted ring model to the 
H I velocity field extracted from the RW H I data cubes (black 
filled circles). Error bars represent the r.m.s. spread of the veloc- 
ities in a single ring of the tilited ring model fitted to the RW 
velocity field. For the sake of clarity, error bars are not shown for 
the model SI rotation curve. 

velocity slices extracted from the Hi data cube. The over- 
lays for tilted ring model SI are presented in Fig 1121 The 
projected rotation curves of the two tilted ring models are 
almost identical to one another with model SI fitting well the 
high-intensity regions of the position-velocity slices. Given a 
set of fixed radial profiles, the ROTCUR routine will adjust 
the rotation curve so that the line-of-sight velocities best 
match the data. Thus although the circular rotation curves 
(i.e. panel F of Fig. llOp differ, the good agreement between 
the projected rotation curves is not surprising. 

Using the smoothed, fixed radial parameter profiles of 
model SI and ring widths of dr — 10", a tilted ring model 
was also fitted to the third-order Gauss-Hermite velocity 
field extracted from the RW Hi cube. The resulting rota- 
tion curve is compared to the rotation curve of model SI 
in Fig. [13] The figure demonstrates that for R < 250", the 
derived rotation curves are very similar to one another. Be- 
yond R ~ 250" the filhng factor of the RW H I velocity field 
becomes too low for meaningful derivations of the rotation 
velocity. For these reasons, only the rotation curves of tilted 
ring models CI and SI are used as mass model inputs in 
SecH 



5 MASS MODELING 

Dwarf galaxies such as NGC 2915 serve as useful probes 
of dark matter (DM). In this section we derive two mass 
models for NGC 2915, one for each of the tilted ring models 
CI and SI. 

The gravitational potential within a galaxy is deter- 
mined by the combined gravitational potentials of all the 
meiss components. The derived 3.6 |xm surface brightness 
profile of NGC 2915 shows little deviation from an exponen- 
tial light distribution at any radius. We therefore concluded 
that the galaxy does not contain a significant central stellar 



bulge. The total mass was treated as the sum of the masses 
of the stellar and gas disks as well as the DM halo. Since 
M oc at a given R, the rotation curves of the individual 
mass components, when summed in quadrature, will yield 
the square of total rotation curve, Vtot'- 

VL = C,,as V,l, -f T, + V^M , (6) 

where Vgas, V*, and Vdm are the contributions to the total 
rotation curve of the gas, the stars and the DM respectively. 
Vtot should match the observed rotation curve. T» is the 
stellar mass-to-light ratio used to convert the observed dis- 
tribution of star light to a stellar mass distribution while 
ctgas ~ 1.37 is the scaling factor used for all mass models to 
take into account the contribution of helium to the rotation 
curve. 



5.1 Gas and stellar distributions 

We used the GIPSY task ROTMOD ()Casertand [l983l ) to 
convert the observed mass distributions into rotation curves. 
Each of the fitted tilted ring models (Sec. I4.3]l was used to- 
gether with the Hi total intensity map to yield an Hi sur- 
face density profile. These profiles were used as ROTMOD 
input which, assuming an infinitely thin gas disk, yielded 
two different Hi rotation curves. Ellipses fitted to 3.6 [im 
surface brightness isophotes suggest a constant inclination 
of i ~ 55° for the old stellar population. Using this inclina- 
tion and a position angle of 306°, a single surface density 
profile for this mass component was derived. This surface 
density profile was extrapolated out to _R = 510" to match 
the radial extent of the observed rotation curve, and then 
converted into a s tellar rotation curve assuming a sech^ ver- 
tical distribution (|van der Kruit|[l98ll ) and h/zo — 5, yield- 
ing zo — 0.12 kpc. 

In principle the DM density distribution of a suitably 
rotating system can be determined knowing the contribu- 
tions to the observed rotation curve from the gas and the 
stars. The Hi rotation curve is directly related to the ob- 
served H I distribution, assuming that the H I is optically 
thin. The situation is not as simple for the stellar compo- 
nent of a galaxy, however, since a stellar mass-to-light ratio 
is needed to link the stellar light and mass distributions. 
Determining the stellar mass-to-light ratio from the rota- 
tion curve alone is problematic. Given an observed rotation 
curve, equally well- fitted mass models ca n be obtained for 
a range of stellar mass-to-light ratios (e.g. Ivan Albada et all 
1985; Swators 1999). This degeneracy therefore leads to un- 
certainties in stellar mass-to-light ratios determined solely 
by r otation curves. 

iBeU fc de Jond (1200 ll ) used spectroscopic and photo- 
metric spiral galaxy evolution models to show that a rea- 
sonable correlation exists between the optical and infrared 
stellar mass-to- light ratios and the colours of integrated stel- 
lar populations. They further showed that the relative trends 
between model stellar mass-to-light ratios and colours are 
fairly insensitive to uncer tainties in stellar population and 
galaxy evolution modeling. lOh et al.l l(2008i') used stellar pop- 
ulation synthesis models (together with a scaled Salpeter 
1955 IMF) to determine a relationship between the stellar 
mass-to-light ratios in the 3.6 \im and K bands (T^'^ and 

respectively): 
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rl-^ = 0.92 X Tf - 0.05. (7) 

From lBell fc de Jond l|200ll ). the relation between the B — V 
colours and the K-haxiA mass-to-light ratio is 



logio (Tf ) = X (B - l^) + a^ 



(8) 



where a^^ and h'' for varying colours and metallicities are 
provided in their Table 4. By combining these two equations 
we can estimate the stellar mass-to-light ratio at 3.6 \im. We 
us e = — .59 a nd = 0.60 for B — V colours as used 
bv lOh et ai] l|2008f ) for IC 2574, the justification being that 
IC 2574 is also a low metallicity, star-forming dwarf galaxy 
with B-V ^ 0.6, much like NGC 2915 B-F ^ 0.3±0.2 for 
NGC 2915 l|de Vaucouleurs et al.|[T99ll : ISchlegel et al.lll998l ') 
thereby yielding (from Eqn.|B]) « 0. 39. Substitut i ng thi s 
value into Eq n. [7| yields » 0.3. iLerov et all (|2008h . 

applying the iBeU et al.l (|2003h relation between B — V 
colour and Tf" and assuming a iKroupal l|200ll ) IMF, found 
Tf = 0.48 - 0.60 for the Jf-band steUar mass-to-light 
ratios of all of the THINGS galaxies in their sample. 



5.2 Dark matter halos 

The way in which dark matter is distributed on galactic 
length-scales has been a topic of much debate. Over the past 
decade, the constant increase in computational power has 
allowed detailed A^-body simulations of the clustering prop- 
erties of cold dark matter particles to be carri ed out. One of 
the fir st large-scale studies was carried out bv lNavarro et al] 
l|l996t ) in which they studied the equilibrium density profiles 
of dark matter halos in a hierarchial galaxy clustering sce- 
nario. Their results suggested the existence of a "universal 
dark matter density profile" l|Navarro et al.|[l997l ) which is 
independent of clustering scale, mass or size as well as the 
power spectrum of initial fluctuations. Their parameterisa- 
tion of this universal profile is 



p{r) 



pcrit 



{r/rs){l + r/r,Y 



(9) 



where 5c is a measure of the density of the Universe at the 
time of collapse of the DM halo, is the characteristic scale 
radius and pcrit ~ 3H^ /SivG is the critical density required 
for closure. This profile scales as r", with a = — 3 and 
Q = — 1 at large and small radii respectively, thereby pre- 
dicting extremely steep inner density profiles of DM halos, 

known as central cusps. 

The Aquarius Project (jSpringel et al.l l2008l ) recently 
completed a set of high-resolution numerical simulations to 
study the formation and structure of galajcy-sized DM h alos 
in the standard ACDM cosmology. 'Navar ro et al.l l|2008l ) re- 
ported that the spherically averaged density profile of the six 
different-sized simulated galaxy halos becomes progressively 
shallower inwards and, at the innermost resolved radius, has 
a logarithmic slope a — dlnp/dlnr > — 1. These predictions 
of cuspy inner DM density profiles can be tested against 
real, dark-matter-dominated galaxies. If cuspy NFW halos 
are found, constraints can be put on halo concentrations and 
cosmological parameters. The rotation curve resulting from 
an NFW density profile is 



Vnfw \ ^ 
, V200 / 



ln(l + cx) — cx/{l + cx) 
xln{l + cx) - c/(l + c) ' 



where V200 is the circular rotation speed at 7-200, the radius 
at which the density of the DM halo equals 200 times the 
critical density, pcrif, x = r/r2oo is the radius in units of 
the virial radius and c — r^oo/rs is the "concentration" pa- 
rameter of the halo. 

Low surface brightness dwarf galaxies are thought 
to be dark-matter-domin ated down to small radii 
(|de Blok fc McGaugh|[l997l ) and can therefore be more di- 
rectly compared to the cold dark matter simulations. When 
this is done, it is often the case that the density profiles of 
simulated DM halos are too steep to fit the rotation curves 
of dwarf galaxies (|N 
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The 

distribution of the DM in these dwarfs is instead consistent 
with approximately constant inner dark matter densities. 
A particular dark matter density profile that describes 
the DM distribution of d warf systems well is that of the 
pseudo-isothermal sphere l|Binnev fc Tremainelll987l ): 



p{r)iso = po 



l+(^ 

Tc 



(11) 



(10) 



where po is the central DM density and rc is the core radius. 
This parameterisation, which has no particular physical jus- 
tification, has a constant density core. The rotation curve 
resulting from such a density profile is 

Viso = (iTrGporc 1 - ^ SiTctan \ ' . (12) 



5.3 Fitted models 

We used the GIPSY task ROTMAS to construct mass mod- 
els for NGC 2915. The task subtracts from the observed 
rotation curve the scaled rotation curves of the stellar and 
gas disks. It then fits to the residuals a rotation curve corre- 
sponding to a particular parameterisation of the DM halo. 
When matching the total and observed rotation curves, the 
points were uniformly weighted to ensure that each point 
in the rotation curve contributed equally to the fit. This 
essentially reduced the fitting procedure to a non-weighted 
least-squares fitting method. 

As previously mentioned, each of the rotation curves 
from the two tilted ring models fitted to the H I velocity field 
were used as mass modeling input. For each of these mass 
models we carried out fits for a pseudo-isothermal sphere 
and an NFW halo. Upper and lower limits for the stellar 
mass-to-light ratio were determined by scaling the contribu- 
tion of the stellar rotation curve to the total rotation curve. 
Under the so-called "maximum disk assumption", the stel- 
lar rotation curve is scaled up to contribute maximally to 
the observed rotation curve at inner radii. Such a scaling 
usually allows the stell ar disk to explain most of the inner 
rotation curve (van Alb ada et al.l fl985: Swatcrs 1999). Un- 
der the "minimum disk assumption" , the contribution of the 
stellar disk is made as low as possible (often zero) while still 
allowing a good fit of the total rotation curve to the observed 
rotation curve. Each of the halo fits was divided into 3 sub- 
cases: A fit 1) with Tt fixed to our predetermined value of 
0.3; 2) using a maximum disk assumption; and 3) using a 
minimum disk assumption. For the minimum disk cases, T« 
was fixed to zero (no contribution from the stellar disk to 
the observed rotation curve). For the maximum disk cases. 
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regardless of the DM halo parameterisation, T* was fixed to 
0.64 and 0.44 for mass models SI and CI respectively. In all 
models, Ugas was fixed to 1.37. Thus, six mass models were 
fitted to each of the 2 rotation curves derived from the tilted 
ring modeling. 



5.4 Results 

The fitted parameters for all of the mass models are sum- 
marised in Table [3] The results for mass model SI are also 
presented in Fig. [I4l This figure shows that NGC 2915 is 
dark-matter-dominated at nearly all radii (depending on 
modeling assumptions). We calculate mass ratios as well as 
the total mass to B-band light ratio at the last measured 
point of each fitted total rotation curve. For mass models SI 
and CI, DM constitutes, on average, 92.8% of the dynami- 
cal mass. This dark-matter-dominance is confirmed by the 
total-mass-to-light ratio of NGC 2915 which is 141 Mq/Lq 
and 129 Mq/Lq on average for mass models SI and CI re- 
spectively. These estimates a re almost double th at of the 
76 Mq /Lq upper hmit set by iMeurer et"al] l| 19961 ) thereby 
making NGC 2915 one of the most dark-matter-dominated 
late-type dwarf galaxies known. 

The second main mass-modeling result that Fig. [14] 
demonstrates is that, for either parameterisation of the DM 
halo, the steeply rising portion of the observed rotation 
curve is poorly matched by the total rotation curve. It is 
usually the stellar rotation cur ve that is used to match the 
inner observed rotation curve (jSwaterj|l999l), yet no such 
scaling is useful in the case of NGC 2915 due to the intrin- 
sically different shapes of the stellar and observed rotation 
curves. The stellar disk is clearly contained well within the 
steeply rising portion of the observed rotation curve which 
continues to rise as V{R) oc R out to J? ~ 150". The best- 
fitting mass models are those in which the stellar rotation 
curve contributes zero to the observed total rotation curve. 
Figure [TJ] shows that when the fitting routine attempts to 
fit the inner observed rotation curve with the DM rotation 
curve, the result is an over-estimation of the observed rota- 
tion velocities for i? < 2 kpc (110") and an under-estimation 
for radii 2 kpc ^ R ^ 5.5 kpc. The poor matches between 
the total and observed inner rotation curves do not allow a 
particular DM halo parameterisation to be confidently ruled 
out or confirmed. Furthermore, the observed non-Gaussian 
line profiles near the centre of NGC 2915 (Sec. I3.6|l lead 
to uncertainties in the shape of the inner portion of the 
observed rotation curve. These, in turn, lead to associated 
uncertainties in the shapes of the inner portions of the fitted 
DM halos. Assuming that the stars are less susceptible than 
the gas to the effects of non-cricular motions, spectroscopic 
observations of the stars would have to be carried out to 
better constrain the inner rotation curve. 

What is the explanation for the discrepancy between 
the observed and total rotation curves? Have we failed to 
include a significant mass component during the decompo- 
sition process? A stellar bulge cannot be the answer since 
our derived 3.6 [im surface brightness profile shows no sig- 
nificant deviation from an exponential light distribution. A 
significant mass of molecular gas at inner radii is a possi- 
bility, yet CO e mission in low-mas s galaxies is usuall y weak 
or not detected (jTavlor et al.ll 19981 : iLerov et al.ll2005l ). How- 
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Figure 15. Mass modeling results for the minimum disk and 
T» = 0.3 cases in which the DM halo is modelled as a pseudo- 
isothermal sphere, c ompared to the mass modeling results of the 
THINGS sample of Ide Blok et al.l l|2008h . The black-filled and 
open triangles correspond to the minimum disk results for mod- 
els SI and CI respectively while the black-filled and open squares 
correspond to the T* = 0.3 results for models SI and CI re- 
spectively. The mass modeling results from THINGS are shown 
as open circles for dwarf galaxies and filled circles for all other 



ever, the conversion factor used to convert CO luminosity 
to molecular gas mass is largely uncertain. 

If we consider the mass-modeling results as they ap- 
pear in Table [3] then the pseudo-isothermal sphere al- 
lows for the best match to the observed rotation curve. 
The fitted pseudo-isothermal sphere, under the minimum 
disk assumption, has a core density and a core radius of 
= 0.17 ± 0.03 Mq pc"^ and = 0.9 ± 0.1 kpc for 



Po 



model SI. These results are similar to those of lMeurer et al.l 
1 19961 ) who, for their favoured mass model, estimated po = 
0.10 ± 0.02 Mq pc"^ and = 1.23 ± 0.15 kpc. We also 
fitted the mass models with T, as a free parameter yet al- 
ways found the best-fitting value to be negative , cons istent 
with the mass modeling results of iMeurer^ et al.1 (|l996l ). 

In Fig. [TS] we compare our derived pseudo-isothermal 
halo parameters for N GC 2915 to t hose derived for the 
THINGS sample from Ide Blok et al.1 (|2008 ). The infrared 
stellar mass-to-light ratios used in their mass models were 
derived from the 3.6 \im images in combination with stellar 
population synthesis arguments. Our models for the mini- 
mum disk and T, = 0.3 cases for both models SI and CI 
are shown in Fig. 1151 This comparison suggests that NGC 
2915 has a central DM core that is very compact and dense 
compared to other late-type systems. 



6 SUMMARY AND CONCLUSIONS 

We have obtained new deep, high-resolution Hi synthesis 
observations of NGC 2915 and have used them to carry a 
detailed study of the H l distribution and dynamics of this 
galaxy. The galaxy has very different optical and Hi prop- 
erties. Optically, it is similar to other dwarf galaxies, con- 
sisting of a small stellar disk. Surrounding the stellar disk, 
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Table 3. Mass modeling results 



1 2 3 4 5 6 

DM halo ISO ISOmax ISOmin NFW NFWmax MFW^in 



Model SI 

1 Xred 38.1 69.8 28.2 62.5 93.4 44.7 

2 r.m.s. (km s~^) 0.66 1.30 0.79 0.75 1.37 0.40 

3 T, (Mq/Lq) 0.3 0.64 0.0 0.3 0.64 0.0 

4 po (10^^ Mq 100.6 ± 23.3 56.9 ± 17.8 178.7 ± 33.0 

5 Rc (kpc) 1.3 ± 0.2 1.7 ± 0.4 0.9 ± 0.1 

6 V2oo(kmB~^) ... ... ... 69.6 ± 7.9 83.5 ± 20.1 62.7 ± 4.1 

7 C ... ... ... 11.7 ± 2.2 7.9 ± 2.7 15.6 ± 2.0 

At last measured pt: 

8 Mtot (10** Mq) 151.4 161.3 161.6 161.1 154.2 148.3 

9 MHi/Mtot (10~^) 43.84 43.86 43.79 43.92 43.03 44.76 

10 Mom /Mtot (10~^) 92.94 89.91 96.62 92.93 90.10 96.62 

11 Mdm/Mj„„ 6.68 4.47 21.84 6.67 4.67 21.34 

12 Mtot /Lb (Mq/Lq) 140.72 140.69 140.91 140.47 143.37 137.9 



1 x;ed 46.7 53.6 37.5 78.1 89.1 61.2 

2 r.m.s. (km s-^) 0.74 1.00 0.27 1.03 1.28 0.66 

3 T, (Mq/Lq) 0.3 0.44 0.0 0.3 0.44 0.0 

4 PO (lO"-'' Mq pc"-') 127.9 ± 34.9 103.2 ± 31.1 227.6 ± 50.1 

5 flc (kpc) 1.10 ± 0.2 1.2 ± 0.2 0.8 ± 0.1 

6 V2oo(kmB~^) ... ... ... 66.4 ± 8.3 70.3 ± 11.5 60.7 ± 4.7 

7 C ... ... ... 12.9 ± 2.8 11.2 ± 2.9 17.0 ± 2.6 

At last measured pt: 

8 Mtot (uf Mq) 138.8 138.3 139.8 137.6 138.4 136.3 

9 Mni/Mtot (10^^) 47.29 47.43 46.92 47.66 47.42 48.16 

10 Mj^M/Mtot (10"^) 92.34 91.13 95.31 92.28 90.92 95.18 

11 M^M/Mi^rn 6.11 6.16 20.31 6.06 5.01 19.77 

12 Mtot/^B (Mq/Lq) 129.00 128.62 130.01 127.99 128.64 126.68 

Comments on rows: Row 1: Reduced goodness-of-fit statistic; Row 2: r.m.s. of the difference between the observed and total 
rotation curves; Row 3: stellar mass-to-light ratio; Row 4: central density for pseudo-isothermal sphere; Row 5: core radius for 

pscudo-isothcrinal sphere, Row 6: circular rotation speed at virial radius for NFW halo; Row 7: concentration parameter for NFW halo; 
Row 8; dynamical mass; Row 9: Hi to total mass ratio; Row 10: Dark to total mass ratio; Row 11: Dark to luminous mass ratio; 
Row 12: Total mass to B-band light ratio. 

Comments on columns: For each model, columns 1 and 4 show the mass modelling results for the case in which T» is fixed to the 
predetermined value of 0.3. Columns 2 and 5 show the results for the maximum disk case while columns 3 and 6 are for the minimum 
disk case. 



however, is a huge Hi disk extending out to galactocentric 
radii greater than R — 510" with wcU-dofinod spiral struc- 
ture in its outer parts. We estimated the total mass of the 
Hi disk to be Mhi ~ 4.4 x 10* Mq. 

Our high resolution, high sensitivity observations 
clearly resolve the inner H I distribution of NGC 2915 into 
two Hi concentrations separated by ~ 1.1 kpc. This obser- 
vation discredits the claim by previous investigators that the 
central gas disk of the system is made up of a massive H I bar. 
Wo have shown the shapes of H I line profiles to bo signifi- 
cantly non-Gaussian in the vicinity of these central H I over- 
densities, being split by ~ 30 km s"'^ on average over the 
spatial extent of the lower Hi concentration. These profiles 
are indicative of non-circular velocity components within the 
gas near the centre of NGC 2915. In our forthcoming paper 
(Elson et al., 2010, in prep.) we investigate in further de- 
tail the central gas dynamics of this system. Linked to these 
non-circular velocity components, we have provided specu- 



lative evidence for the possibility of cold gas accretion from 
the inter-galactic medium onto the outer disk of NGC 2915 
which, if confirmed, would have significant implications for 
the evolutionary history of this system. These complicated 
dynamics lead to the now picture that NGC 2915 is not a 
simple, isolated, low-surface-brightness galaxy but rather an 
evolving system with a complex interplay between its vari- 
ous mass components. 

The main focus of this paper has been the regular H I 
dynamics of NGC 2915. We have fitted a tilted ring model 
to the Hi velocity field. The best-fitting model was one in 
which the H I disk is severely warped, with the inner and 
outer disks inclined at i ~ 70 - 75 ° and i ~ 55° respectively. 
This result is, however, made uncertain by the non-Gaussian 
line profiles near the centre of the galaxy. We therefore also 
fitted a model with an almost constant inclination of i ~ 55° 
for the entire H I disk. Both models yielded a rotation curve 
typical of a late-type spiral with V{R) oc R rotation out to 
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Figure 14. Pseudo-isothermal sphere (ISO) and NFW rotation curve fits for NGC 2915. The panels on the left show the fits using the 
pseudo-isothermal sphere while the panels on the right show the fits using an NFW halo. The top panels show the fits when the stellar 
rotation curve is scaled by our predetermined value of T« = 0.3 while the middle and bottom panels shows the fits for the maximum 
and minimum disk cases respectively. In all panels the black filled circles represent the observed rotation curve while the error bars 
represent the r.m.s. spread of velocities in a given ring. The black dot-dash curve shows the rotation curve of the gas while the black 
dotted curve shows the rotation curve of the stellar disk. The solid grey curve shows the resulting rotation curve of the DM halo. Finally, 
the solid black curve represents the total rotation curve resulting from the best-fit model. In each panel the goodness-of-fit statistic is 
presented. The parameter radial profiles used to constructed the stellar and gas rotation curves were those of tilted ring model SI. The 
observed rotation curve is also from tilted ring model SI. 
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R ~ 150" and constant velocity thereafter. Rotation curves 
separately derived for the approaching and receding sides of 
the galaxy differed significantly within their steeply rising 
portions, providing the clear evidence that NGC 2915 is a 
kinematically lopsided galaxy. 

Our results show that these rotation curves cannot be 
accounted for by the stellar or gas disks of NGC 2915. By 
using the observed rotation curves as mass model inputs, we 
found the galaxy to be dark-matter-dominated at nearly all 
radii with the stellar disk unable to account for the V(R) oc 
K portion of the rotation curve. The best-fitting mass model 
is one in which the stellar disk does not contribute to the 
observed rotation curve and in which the inferred dark mat- 
ter halo is parameterised as a pseudo- isothermal sphere. This 
model has a central core density of po = 0.17±0.03 Mq pc~'^ 
and a core radius of = 0.9 ±0.1 kpc. The non-circular 
gas motions at inner radii lead to uncertainties in the fitted 
dark matter halo parameters. These, together with the fact 
that neither DM halo parameterisation allows for an accu- 
rate match of the observed rotation curve at inner radii, do 
not allow a particular DM halo parameterisation to be con- 
fidently ruled out or confirmed. Further observations and 
modeling are required to more accurately determine the in- 
ner rotation curve of NGC 2915 and to therefore more con- 



fidently discriminate between various parameterisations of 
the DM halo. It is clear, however, that NGC 2915 has a very 
dense and compact dark matter core. At the last measured 
point on the rotation curve, this galaxy has a total-mass-to- 
light ratio of MjL ~ 140 Mq/Lb, almost twice that of the 
upper limit placed mv lMeurer et al] l| 19961 ). thereby making 
it one of the most dark-matter-dominated galaxies known. 
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